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ABSTRACT 



An analysis of steady state operation of a vertically 
oriented, variable conductance heat pipe is presented. The 
effects of binary mass diffusion, axial pipe wall conduction, 
and gravitational effects, caused by the difference in 
molecular weight of the non-condensible gas and the working 
fluid are incorporated. Analytical expressions for the con- 
servation of mass, momentum, and energy are combined along 
with equations of state to describe steady state operation 
of the heat pipe. These expressions are combined to form 
a system of three differential equations with three unknowns; 
working fluid vapor velocity, mass fraction, and temperature. 
These equations are nondimensionalized and an iteration scheme 
for numerical solution of the equations is presented . 
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INTRODUCTION 



I . 



A. BACKGROUND 

A heat pipe is a closed, evacuated chamber that has the 
inside walls lined with a wick structure, saturated with a 
volatile working fluid. The steady state operation involves 
vapor heat transfer and capillary action. A comprehensive 
coverage of the theory and operation of heat pipes can be 
found in Dunn and Reay [Ref. 1], Briefly, the operation of 
a vertically oriented heat pipe will be discussed. The heat 
pipe has three sections: an evaporator section, an adiabatic 

section, and a condenser section. The heat pipe with the 
three sections marked is shown in Figure 1.1. The lower evapora- 
tor end of the heat pipe is heated, causing the working liquid 
in the wick to vaporize. This phase change causes the working 
fluid to absorb latent heat and the vapor pressure increases 
in the evaporator end. The working fluid is drawn away from 
the wick, causing the liquid pressure to decrease. The vapor 
moves to the cooler condenser where it is condensed and the 
working fluid gives up its latent heat. The condensate is 
deposited on the wick which causes the liquid level to in- 
crease thus increasing the liquid pressure. Capillary forces 
take the higher pressure liquid from the condenser to the lower 
pressure liquid end of the evaporator. The capillary forces 
are assisted by gravitational forces when the heat pipe is 



9 



Q 



out 



Q in 




x = 0 



Condenser 



▼ 



x = L 



Adiabatic Section 



Evaporator 



Figure 1.1. Heat Pipe 
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vertically oriented as shown in Figure 1.1. it is important to 
note that for heat pipes the evaporator position is not re- 
stricted, but since the experimental data and numerical model 
use a vertically oriented heat pipe, this orientation will be 
used throughout this thesis. 

The great advantage of heat pipes is that the vaporization 
and condensation process allow for high heat transfer rates 
with small temperature differences. Two other advantages are 
that the basic heat pipe has no moving parts and can be used 
in many geometrical configurations. The heat pipe has limi- 
tations and again, a comprehensive presentation of these limi- 
tations can be found in Dunn and Reay [Ref. 1] . Briefly, the 
five heat pipe limitations are as follows : 

(a) The viscous limit is basically the low temperature limit. 
This would be the minimum temperature for the heat 

pipe to begin operation. Generally, the viscous limit 
is only important during start up of some heat pipes, 
i.e., liquid metal pipes. 

(b) The sonic limit describes the limit of the vapor flow 
velocity at the heat pipe evaporator exit. The re- 
sulting choked flow condition limits the axial heat 
flux capacity of the heat pipe. 

(c) The entrainment limitation is due to the removing of 
the liquid before it can return to the evaporator by 
vapor flowing at high velocity. 

(d) The wicking or capillary limit reflects the fact that 
each combination of wick and fluid has a maximum 
capillary head. As stated earlier, the capillary head 
must increase with the liquid and vapor pressure drops 
and can do so only up to the maximum capillary head. 

(e) The last limit is the boiling limit. If temperatures 
in the evaporator get too high the working fluid will 
begin to boil within the wick. This will disrupt liquid 
flow to the evaporator, causing the heat pipe to stop 
working . 
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These five heat pipe limitations vary in importance, depending 
on the temperature and corresponding axial heat flux operating 
point. Since at different temperatures one can encounter 
different heat pipe controlling limitations, it is important 
to have a basic understanding of all five limitations. 

B. VARIABLE CONDUCTANCE HEAT PIPE 
1 . Basic Operation 

Variable conductance heat pipes are described in detail 
in Marcus [Ref. 2]. As described in the background section, a 
change in the heat load will automatically change the operating 
temperature and pressure of the heat pipe, given a constant sink 
temperature. This is an undesirable quality for many situations 
where it is desired to have a stable source temperature under 
conditions of varying heat load or sink temperatures. An 
example of this case would be the cooling of sensitive elec- 
tronic equipment. Basically, a noncondensible gas is introduced 
into the heat pipe along with the working fluid, making the 
heat pipe become a passively controlled device. The nonconden- 
sible gas occupies a certain portion of the condenser section, 
depending on operating conditions. Increasing operating tempera- 
ture in turn increases the vapor pressure of the working fluid 
which compresses the noncondensible gas into a smaller volume. 
This provides a greater active condenser area, thus allowing 
a more stable source temperature than a conventional heat pipe. 
When the operating temperature decreases, working fluid vapor 
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pressure falls and the noncondensible gas expands, blocking a 
larger portion of the condenser. The result is a passively 
controlled variable condenser area, which in turn makes a 
more stable source temperature. 

2 . Variable Conductance Heat Pipe Theory 

Marcus [Ref. 2] describes a flat front theory and a 
'diffuse front theory' for variable conductance heat pipes. 

Flat front theory assumes the noncondensible gas is in the 
form of a plug with sharply defined vapor-gas interface which 
is perpendicular to the axis of the pipe. Since this model 
neglected the axial conduction in the heat pipe wall and diffu- 
sion of vapor into the stagnant gas plug, experimental results 
[Ref. 2] have shown this model to be inaccurate. The early 
model of flat front theory provided initial insight into the 
operating characteristics of gas loaded heat pipes. With the 
incorporation of binary diffusion between the working fluid 
vapor and the noncondensible gas and axial conduction in the 
pipe wall, Marcus' [Ref. 2] 'diffuse front theory' establishes 
a model of gas loaded heat pipes which provides better simi- 
larity between computer program and laboratory results. 

Gravitational effects on the vapor-gas interface are not 
included in flat front theory or the 'diffuse front theory.' 
Gravitational effects can greatly distort the axial temperature 
profile when the working fluid and noncondensible gas have 
significantly different molecular weights. This fact is demon- 
strated by experimental results reported by Kelleher [Ref. 3]. 
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Batts [Ref. 4] conducted variable conductance heat pipe experi- 
ments using working fluid/noncondensible gas combinations of 
methanol/helium and methanol/krypton . Kane [Ref. 5] developed 
an analytical model of variable conductance heat pipes which 
includes the effects of binary mass diffusion, axial pipe-wall 
conduction, and gravitational effects which exist at the vapor- 
gas interface region as shown by Batts [Ref. 4] experiments. 
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II. 



FORMULATION OF AN ANALYTICAL MODEL 



A. PROBLEM STATEMENT 

The purpose of this chapter is to formulate the equations 
describing the steady state operation of a vertically oriented, 
gas loaded, variable conductance heat pipe. As in Figure 1.1, 
the condenser is above the evaporator for this orientation. 

This analysis includes the effects of binary mass diffusion, 
axial pipe wall conduction, and gravitational effects due to 
the difference in molecular weight between the noncondensible 
gas and the working fluid. 



B. ASSUMPTIONS 

The following assumptions are used to develop the system 
of equations: 

1. Heat, mass and momentum transport processes are treated 
as one-dimensional in the axial direction within the 
vapor space of a vertically oriented heat pipe. 

Kelleher [Ref. 3] shows that in any other orientation, 
gravitational effects cause highly three-dimensional 
spatial variations when the molecular weights of the 
working fluid and noncondensible gas differ substantially. 

2. The noncondensible gas acts as a perfect gas. The 
mixture of noncondensible gas and working fluid are 
treated as thermodynamically independent substances. 

3. Because of direct contact between the vapor and the 
liquid saturated wick, the working fluid vapor is at 
saturation conditions at all times. 

4. The thermal resistance of the wick-wall structure, in 
the radial direction, is negligible in comparison with 
the thermal resistance due to convection outside the 
pipe. 
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C. ANALYSIS 



Kane [Ref. 5] presents a proposed solution for this prob- 
lem. The following is a significant change to his proposed 
solution with modifications being noted as they occur in 



this section. Variables with the subscript 'a 1 


' apply to the 


working fluid and variables with the subscript 
the noncondensible gas. 


'b' apply to 



Utilizing the assumption that the noncondensible gas is 



a perfect gas, the equation for noncondensible 
is : 


gas pressure 


P = o — T 

b M b M fa 


(II. 1) 



The principles of mixtures of thermodynamically independent 
substances, which are found in Reynolds and Perkins [Ref. 6], 
are used to define the total density and the total pressure. 



D = 0 + 0, 

H y a y o 


(II. 2) 


P = P + P, 
a b 


(II. 3) 


The definition of mass fraction is used to 
for the working fluid vapor mass fraction. 


form the equation 


_ ^a 
W a “ p 


(II. 4) 
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The solution, utilizing the conservation of mass, is 
based on a differential element representing a section of 
the heat pipe condenser, Figure 2.1. The conservation of 
mass of species 'a' for the differential element yields: 



V 

a a c 



(p V + ^— (p V )dx)A 
M a a dx 'a a' c 



+ m P dx . 
a v 



This equation sets the mass of species 'a 1 into the control 
volume equal to the mass of species 'a' out, including a 
term to account for the removal of species 'a' by condensation. 
The process of condensation is accounted for by multiplying 
the condensation mass flux times the area, defined by the 
vapor space perimeter times the length, dx, of the differen- 
tial element. After basic algebra, the equation for conser- 
vation of mass of species 'a' becomes: 



0 



m P 
a v 



A ^-vp V ) 
c dx a a 



By using Fick ' s law of binary diffusion [Ref. 7], and solving 
for p V , the resulting formula is: 

cl a 



P V 

a 



dm 

“a oV - pD -a# 



The final form for the conservation of mass for species 'a' 
becomes : 



. P . , dw 

v , d , lT . d , _ a 

m a + dJ ( “a pV) - d? (pD 



dx 



) = 0 



(II. 5) 
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Figure 2 . 1 . Differential Element for Conservation of Mass 



Since species 'b' is a noncondensible gas, the condensation 
mass flux of species ' b ' equals zero. The resulting equation 
for the conservation of mass for species ' b ' is: 





0 . 



The mass continuity of the mixture is derived by adding 
the conservation of mass of species ’a’ with the conservation 
of mass of species ' b ' . 



the final result of the mass continuity of the mixture is: 



Figure 2.2 is a differential element of the vapor space in 
the condenser and is used to derive the conservation of linear 
momentum equation for the mixture of noncondensible gas and 
working fluid vapor. The conservation of linear momentum for 
the mixture of noncondensible gas and working fluid vapor is 
derived by setting the summation of external forces acting on 
the control volume in the x direction equal to the net efflux 
of linear momentum. 




0 . 



Because (o + oo, is equal to one, and (to +co, ) = (1) = 0 

a b ^ dx a b dx 




0 . 



(II. 6) 



I? = / V(pV dA) 

CS,x 
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Figure 2.2. Differential Element for Conservation of Momentum 



The summation of external forces acting in the x direction are 




p g A dx - x P dx - dp A 
c v c 



These terms represent the force due to gravity acting on the 
control volume minus the forces due to wall shear stress and 
pressure acting on the control volume. The wall shear stress 
can be calculated by: 



x 




V 



2 



where A is a dimensionless friction factor. If developed 

1 6 

laminar flow is assumed the Darcy friction factor, A = — , 
could be used with the Reynolds number defined as 



Re 



p VD 
_ v 

y 



The equation for wall shear stress then becomes: 



T 




V 



where y is the viscosity of methanol vapor, V is the velocity 
of the vapor and is the diameter of the vapor space. The 
terms to represent the net efflux of linear momentum are: 



/ V(pvdA) 

CS ,x 



* 2 

V m P dx + (pv +d(pv) (v+dv))A - pV A 

3. 3. V CO 

V m P dx + A pvdv + A vd(pv) . 

a a v c K c 
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By setting the summation of external forces in the x direction 
equal to the net efflux of linear momentum and combining like 
terms, the result is: 



dP 

dx 




(II. 7) 



The above equation represents conservation of linear momentum 
for the mixture of noncondensible gas and working fluid vapor. 
Kane [Ref. 5] assumes zero velocity for the noncondensible gas 
and derives an equation for conservation of linear momentum 
of species 'b.' The species 'b' conservation of linear momen- 
tum equation is not necessary to derive the model in this 
analysis. Since the conservation of linear momentum for the 
mixture of the noncondensible gas and working fluid vapor is 
the more consistent approach. Equation (II. 7) is used through- 
out this thesis. 

Figure 2.3 is a differential element of the condenser wall 
and represents the wick and pipe wall combination as one ele- 
ment. The wick-wall combination is treated like a fin and is 
used to derive the equation for the conservation of energy. 
Heat transfer is represented by the latent heat deposited on 
the 'fin' by condensation of the working fluid vapor plus the 
heat transfer from Fourier's law of conduction entering the 



' fin. ' 



m 




P dx + Q 



c 
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m 



a 





Figure 2.3. Differential Element for Energy Balance 
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Heat transfer out is represented by Fourier's law of conduction 
leaving the 'fin' plus the energy which is convected away from 
the fin. 



Q 



out 



Q + dQ + h P (T -T )dx 

C C S °° 



The equation for conservation of energy is derived by setting 
heat transfer in equal to heat transfer out of the 'fin.' 



d 2 T 

dx 2 



hP h P 

r^-(T -T ) + v m 

kA °° kA a 



= 0 



(II. 8) 



The experimental data from Batts [Ref. 4] is used to 
correlate this analysis. The working fluid is methanol, 

CH^OH. The working fluid vapor in the heat pipe is in direct 
contact with the working fluid liquid in the wick and is 
therefore considered to be a saturated vapor at the local 
temperature of the wick-wall combination. Eubank [Ref. 8] 
tabulates data for the saturation properties of methanol. 

These tables are used to generate functions of local temperature 
for partial pressure, partial density and heat of vaporization 
of saturated methanol vapor. 

(a) The partial pressure of saturated methanol vapor is a 
function. Fa, of local temperature. 



P = Fa (T ) 
3 . 



(II. 9) 



24 



(b) The partial density of saturated working fluid vapor 



is a function, Sa, of local temperature. 



P 



a 



Sa (T) 



(II .10) 



(c) The working fluid heat of vaporization is a function. 



In summary, the eleven equations which are stated in this 
chapter form the foundation for the analysis of a vertically 
oriented, variable conductance heat pipe. The remaining 
chapters derive a numerical model and iteration scheme based 
on the eleven equations from this chapter. 




h 



fg 




(ii . id 
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Ill . 



NUMERICAL MODEL FORMULATION 



A. VARIABLE ELIMINATION 

The system of eleven equations can ultimately be reduced 

to a system of three equations by eliminating variables. 

Equation (II. 9) is eliminated by using the equation of state, 

Fa(T), for P . The partial density of the saturated working 

fluid vapor, p (Equation II. 10) is eliminated by using Sa(T). 

Hj. (T) , Equation (II. 11), is used to eliminate h c . The non- 
fg fg 

condensible gas pressure, , is eliminated by using the perfect 
gas equation of state (Equation II. 1). The treatment of working 
fluid vapor and noncondensible gas as thermodynamically inde- 
pendent substances gives total pressure as the sum of partial 
pressures and total density as the sum of partial densities. 
Total pressure is eliminated by using the sum of 'partial pres- 
sures, Equation (II. 3). The noncondensible gas density is 
eliminated by solving Equation (II. 2) for and replacing p^ 
with p -p . The definition for mass fraction of the working 

3 

fluid vapor, Equation (II. 4), is used to replace total density. 

After making the preceding seven substitutions, four equa- 
tions with four unknowns are left. The four remaining depen- 
dent variables which are functions of axial position are: 

• _ 2 

(a) m , working fluid condensation flux [lbm/ft sec] 

3 

(b) V, velocity of the mixture [ft/sec] 

(c) co , working fluid mass fraction [0 < co < 1] 

(d) T, temperature [°R]. 
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The four remaining equations are: 



(a) m v) 

a A dx to ' 
c a 



= 0 



(III.l) 



b m — ^ + i-(Sa V 
a A dx 
c 



d , SaD dc °a . 

dx (i) dx 
a 



= 0 



(III. 2) 



(c) + ^-(Sa(— -1) -5- T) 

dx dx w M, 

a d 



O n 0 PmP Q 

/Sa x d ,Sa TT 2. v T7 a v 8y „ 

<— 9> - d^ ( — V > - iT V ~ — D- V 

a a c a c v 



(III .3) 



(d) 



d 2 T 

dx 2 



P HP. 

h (T -T ) + — m = 0 

KA °°' KA a 

P P 



(III. 4) 



Furthermore, Equation (III.l) can be eliminated by solving this 

equation for working fluid condensation flux, m , and substi- 

a 

tuting that into the three remaining equations. 

The resulting three equations with dependent variables of 
velocity, working fluid mass fraction, and temperature are: 



d(i) 

(a) ( — V) + (Sa V) - -=— ( — D-j — ) = 0 

dx a) dx dx a) dx 

a a 



(III. 5) 



b Sa — -1 t-5- T 

dx dx 'co M, 

a b 



= (— g) - s-(— v 2 ) +— v) 

a) ^ dx a) u) dx u 

a a a a 



P 8y 
v 

A D 
c v 



V (III. 6) 
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(c) 



0 



(III. 7) 



d 2 T 

dx 2 



KA 



h (T -T ) - 



H A 
fg c 

KA 



V) = 

dx a) 

a 



This set of three equations and three unknowns describe the 
operation of a vertically oriented, gas loaded, variable conduc- 
tance heat pipe. 

B. BOUNDARY CONDITIONS 

These equations describe the thermodynamics of the heat 
pipe condenser. The independent variable, the distance along 
the heat pipe, is measured with the origin at the condenser 
end farthest from the evaporator. The total condenser length 
is defined as L. Based on the system of three equations, two 
boundary conditions are required for temperature and mass frac- 
tion and one boundary condition is required for velocity. 

From experimental data, the temperature at the far end of 
the condenser and in the evaporator are known. The boundary 
conditions for the dependent variable of temperature are: 

when x = 0, T = T , where T equals the temperature at the 

far end of the condenser. 

when x = L, T = T , where T equals the temperature at 

the end of the condenser closest to 
the evaporator. 

The highest velocity of the working fluid vapor is reached 
in the evaporator exit. This velocity can be calculated from 
the known heat input to the heat pipe. The vapor velocity in 
the evaporator, V e , is equal to the heat transfer rate into 
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the evaporator divided by the product of local density of the 
working fluid vapor, cross sectional area of the vapor space, 
and the local heat of vaporization of the working fluid. Due 
to the physical boundary at the end of the pipe, where x = 0, 
the velocity at the farthest end of the condenser away from the 
evaporator is equal to zero. The boundary conditions for the 
dependent variable of velocity are: 



when x = 0 , V = 0 . 0 and 



when x = L, V = V 



Q e 

p a A c h fg 
e J e 



Based on the system of three equations, only one of the above 
boundary conditions is required. The boundary condition at the 
origin is used for the iteration scheme and V g is used to 
normalize velocity. 

The temperature and pressure at the farthest end of the 
condenser from the evaporator are known from laboratory data. 
Assuming thermodynamically independent substances, the total 
pressure is the sum of partial pressures and the total density 
is the sum of partial densities. The saturation pressure of 
the working fluid vapor, P , could be determined from tabulated 
data by entering the saturation pressure table with condenser 
temperature. An example of this table can be found in Eubank 
[Ref. 8] for methanol saturated vapor. The noncondensible gas 
pressure, , can be calculated by subtracting the working 
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fluid vapor pressure from the total pressure. With this data 
the initial mass fractions can be established by using the 
following procedure from Reynolds and Perkins [Ref. 6] : 



X 



a 



P 

a 



P total 



total 



X 



M 



a a 



Lbm 

Lbmole mixture 



for species 'a' 



Xb W b 



Lbm 



Lbmole mixture 



for species 'b' 



0 ) 



X M 

a a 



a 0 x a H a + *b M b 



00 



*b M b 



b. Y M +Xi_Mu 

0 A a a A b b 



At the end of the condenser closest to the evaporator, it is 
assumed that the concentration of noncondensible gas is zero 
so that at this point is equal to one. The boundary condi 
tions for the dependent variable of mass fraction for the 
working fluid vapor are: 

when x = 0 , w = to , and 
a a Q 

when x = L, w = 1.0 . 
a 
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C. EQUATION NONDIMENSIONALI ZATION 

To generalize the equations and facilitate their solution 
it is convenient to use dimensionless variables defined as 
follows . 

(a) The independent variable, condenser length, is defined 
by : 

x 



where x is the distance along the heat pipe, measured from 
the condenser end farthest from the evaporator, and L is the 
total condenser length. 

(b) A dimensionless temperature can be defined by: 




where T is the local temperature, T^ is the evaporator tempera- 
ture, and T^ is the ambient temperature. In this analysis, the 
first and second derivative of temperature with respect to dis- 
tance are required. 



dT _ dT d_0_ d_n_ _ ***e *^°°d6 

dx d0 dr] dx L dri 



and 



d 2 T 

dx 2 



d_,dT. _ d . ^e d9 . drj 

dx dx - dn L dr) dx 



T -T 

e ° 



d^e 

dn 2 
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(c) The second dependent variable, working fluid mass 
fraction, is dimensionless by definition and its numerical 
value is greater than or equal to zero and less than or equal 
to one. The first derivative of working fluid mass fraction 
with respect to distance is: 



doj 

a 

dn 




dn 

dx 



1 

L 



du) 

a 

dg 



(d) Velocity of the working fluid vapor is made dimension- 
less by normalizing local vapor velocity with the evaporator 
vapor velocity, V . The nondimensional velocity and the first 
derivative of velocity with respect to distance are then 
defined by: 



e 

dV _ dV dV* djn = dV* 

dx dV* dn dx L dn 

(e) In order to complete nondimensionalizing process, the 
following other dimensionless variables are used. 

(i) The working fluid vapor density is normalized with 
respect to the evaporator vapor density. 



Sa* 
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(ii) The working fluid vapor pressure is compared 



to the evaporator vapor pressure. 



Fa* 




(iii) The working fluid heat of vaporization is normalized 
with respect to the condenser working fluid vapor heat of 
condensation . 




H 

H 



fg 

fg,c 



(iv) The working fluid vapor viscosity is normalized 
with respect to the evaporator vapor viscosity. 



y * 




After substitution of nondimensional terms, the three 
equations which describe operation of a vertical variable 
conductance heat pipe become: 



(a) - ^-( 



d ,Sa*.V v 



D ,d Sa’ 



dco 



dn ' w > + dn (Sa * y * ) LV ( dn w “ dn 
a e a ' 



) = 0 (III. 8) 



The dimensionless coefficient from Equation (III. 8) is 
D 



N1 = 



LV * 
e 
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Fa 



(b) 



* R(T -T ) , , 

e dFa* , e 00 ,d . , 1 . . . . •» 



Sa V 2 dn 
e e 



^ (sa *^- 1)(9t vf 

b e 



■) ) 



Lg Sa* 

,,2 00 

V a 
e 



d_ ( Sa*V^_ } + v * d^ ( Sa*V* } 



dr| cj. 



dn 03. 



P 8 p L 
v e * 

u* V* 

A D Sa V M 
c v e e 



(III. 9) 



The dimensionless coefficients from Ecruation (III. 9) are 



Fa 



N2 = 



Sa V 2 
e e 



N3 - 



R (T -T ) 

e 00 

m,v 2 

b e 



N4 = 



T -T ' 

e 00 



N5 = 



and N6 
V 



p 8 n l 

v e 

A D Sa V 
c v e e 



,2. p hL 

d 0 s 0 

( } , 2 K A 0 

dn p 



A Sa V L 
fg ,c c e e 



H] 



K A (T -T ) fg 
p e 00 



x d = 0 

dn a) 



(III. 10) 



The two nondimensional terms of interest from Equation (III. 10) 
are 



2 

P hL 

N7 = -§— r and N8 

X\ n 

P 



fg,c 

K A 

P 





V 

e 



L 



Numerical evaluation of the dimensionless coefficients is 
required to determine the relative significance of each term. 
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Two experimental runs from Batts [Ref. 4] are used to evalu- 
ate the dimensionless coefficients. A low powered run 
(Q = 22.9 watts) with methanol as the working fluid, and krypton 
as the noncondensible gas, is chosen. This run demonstrates 
an example of the heavier molecular weight noncondensible gas 
diffusing through a large portion of the condenser. A second 
higher powered run (Q = 70 watts) with methanol as the working 
fluid, and helium as the noncondensible gas, is chosen because 
it demonstrates an example of the lighter molecular weight 
noncondensible gas occupying a small portion of the condenser. 

The convection heat transfer coefficient for the heat pipe is 

Btu 

set equal to 1.0 ~ , representing free convection. A 

hr . ft. °R 

Btu 

system thermal conductivity of approximately 4 . 56 ~ 

hr. ft. °R 

results from a parallel circuit combination of conductive 
thermal resistances for the stainless steel pipe wall and 
methanol filled stainless steel wire mesh wick. Kane [Ref. 5] 

calculated representative binary mass diffusion coefficients 

2 2 
of 0.4 ft /hr for krypton-methanol and 1.9 ft /hr for helium- 

methanol. These values for binary mass diffusion are also 

used in this analysis. The experimental values that are used 

to calculate the dimensionless constants are listed in Table 

I. The values of the dimensionless constants for the two 

experimental runs are listed in Table II. 

Based on the tabulated numerical results. Equation (III. 9) 

is modified to: 
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TABLE I 



NUMERICAL DATA 



Constants for both runs: 

R = 1545.33 ft Lbf/Lbmole °R 

g = 32.2 ft/sec 2 

L = 3.5 ft 

A = .0199 ft 2 
c 

A = .001917 ft 2 
P 

D = .15917 ft 

V 

P = .524 ft 
s 

P = .500 ft 
v 

h = 1.0 Btu/hr ft 2 ° R 

K = 4.5585 Btu/hr ft°R 

Low powered krypton/methanol 

Q = 22.9 watts = 78.1348 Btu/hr 

T = 592 °R 
e 

T = 543 °R 
c 

T = 536 0 R 

CO 

D = 0.4 f t 2 /hr 

= 83.8 Lbm/Lbmole 



Gas constant 

Acceleration of gravity 

Length of condenser 

Cross sectional area of 
vapor space 

Axial cross sectional area 
of wick-wall combination 

Diameter of vapor space 

Periphery of outside of 
pipe 

Periphery of vapor space 

Convection heat transfer 
coefficient 

Thermal conductivity 
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TABLE I (CONTINUED) 



V = .041339 ft/sec 
e 

Fa e = 10.24 Lbf/in 2 

Hfg c = 500.21 Btu/Lbm 

Sa e = 0.054871 Lbm/ft 3 

y e = 224.502 xio" 9 Lbf sec/ft 2 

High powered helium/methanol 

Q = 70 watts = 238.84 Btu/hr 

T = 64 0 0 R 
e 

T = 54 0 0 R 
c 

T = 536°R 

OO 

D = 1.9 ft 2 /hr 

= 4.003 Lbm/Lbmole 

V = .050485 ft/ sec 
e 

Fa e = 28.52 Lbf/in 2 

H. =501.11 Btu/Lbm 
fg,c 

Sa e = 0.14446 Lbm/ft 3 

y e = 243.993 x 10 _9 Lbf sec/ft 2 
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TABLE II 



Low 



Nl 



N2 



N3 



N4 



N5 



N6 



N7 



N8 



DIMENSIONLESS CONSTANTS 



Powered Krypton/Methanol High Powered Helium/Methanol 



LV 



D = 7.679 x 10 4 



Nl = =^- = 2.9869 x 10 3 

V 



Fa 



Sa V 2 
e e 



= 506.355 x 10 



N2 - 



Fa e 

Sa V 2 
e e 



= 359.166 x 10 



R (T -T ) 

e °° 

M, V 2 
b e 



= 19 .4581 x 10 N3 = 



R (T -T ) 

e 00 

M, V 2 
b e 



= 507 . 224 x 10' 



T -T 
e 0 



= 9.571 



N4 = 



T -T 
e 0 



= 5.154 



^ = 65948.4 

V 

0 



N5 = ^ 

V 

e 



= 44218.0 



P 8 y L 
v e 

A D Sa V 
c v e e 



= 14.086 N6 = 



P 8 U L 
v e 

A D Sa V 
c v e e 



= 4.76143 



P hL 
s 

K A 



= 734.553 



N7 = 



P hL 
s 

K A 



= 734.553 



A Sa V L 
f g , c c e e 

K A (T -T ) 
p e 00 



N8 = 



H, A Sa V L 
fg,c c e e 

K A p (T e -Tj 



= 581.362 



= 1008.30 
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Fa 



* R(T -T ) 
e dFa* . e °° 



Sa V 2 
e e 



dn 



M, V 
b e 



2— [ ^T (Sa * ( ^r _1) (e + 



T -T 

e < 



-) ] 



_ Lg Sa 1 



V 



2 w 



(III .11) 



The viscous term and the two inertia terms are neglected be- 

1 3 6 

cause their order of magnitude is 10 compared to 10 and 10° 

for the remaining terms. The nondimensional equations (III. 8), 

(III. 10), and (III. 11) are used to develop the analytical model 

of the vertically oriented, variable conductance heat pipe. 



D. FORMULATION OF THE NUMERICAL MODEL 

Equation (III. 8) is integrated and is solved for the first 
derivative of the working fluid mass fraction. Integrating 
both sides of Equation (III. 8) gives: 



Sa * V * + Sa *V* - HI Sa * 



to 



to 



dco 

a 

dn 



= c 



(III. 12) 



The lower endpoint is used (n = 0) to evaluate the constant 
of integration. When n = 0 > V* at the end of the pipe equals 
zero. The constant of integration then becomes: 



C 



dw 



-N1 



Sa 



dn 



(III .13) 



By replacing the value for the constant of integration and 
solving for dw /dn , the final version of Equation (III. 12) is: 

cl 
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(III. 14) 




-V* 

Nl 



+ 



V*co. 

Nl 



CO 



dco 



Sa* Sa dn 



The derivative in Equation (III. 11) is expanded to get the 
following . 



|^(S a *(^-l)(6 + N4>) 



dSa* 

dn 



(rr- - D (0 + N4) 

00 

a 



dco 

+ Sa* -4 -3-^(0 +N4) 



-_1 

2 

00 



dn 



+ 



Sa* (— -1) 

co 

a 



d0 

dn 



(III. 15) 



The chain rule is applied 



to 



dSa* 
dn ' 



dSa* _ dSa* dSa dT d6 
dn ~ dSa dT d0 dn 



T -T 

e °°.dSa d0- 

Sa ' dT dn 
e 



(III. 16) 



Applying the same procedure to 



dFa* 

dn ' 



the result is: 



dFa* 

dn 



T -T 

e dFa d0 

Fa e ’ dT dn * 



(III. 17) 



By combining terms and solving for the final form for 

equation (III. 11) is: 
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d0 

(a) The coefficient of - 3 — is set eaual to G (w ,9). 

dn “ a 



N2 (T -T ) N 3 (T -T ) , 

„ e 00 dFa , e 00 dSa ,1 x x 

G = Fa dr + si 3t - 1 ’ <e +N4) 



+ N3 (Sa*) (— -1) . 



co 



(III. 18) 



d 0 

(b) Solving for the result is 



d0 r N5Sa* . N3Sa* d a,. . 

- 3 — = l + * — —5 — (9 +N4)]/G(co ,9 

dn co 2 dn a 

a co 



(III. 19) 



Equation (III. 10) is integrated and the result is: 



d9 xn r Q j 
dF ' N7 / edn 



N 8 H* 

fg 



Sa*V* 

03 

a 



C . 



(III. 20) 



Using the boundary condition of n = 0 where V* = 0.0, the 
constant of integration is: 



C = 




(III . 21) 



The final version of Equation (III. 10) which is solved for V* 
is : 



V* 



(^03 a -N7(/8dn)03 a 



d9 
co j — 
a dn 



)/(N 8 H* Sa*) 
0 £g 



(III. 22) 
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The final system of three equations with nondimensional 



variables is: 



dw 



dw 



(a) Equation (III. 14), _iL(v*,oo 

dn a dn 



) : 



dw 
a 

dn 



V*U) 



-V* 

N1 + Nl 



a + w a p 0 dw a 



S * Sa dn 
Q e 



fl 

(b) Equation (III. 19), ^( w a , 0 , -g^-) ; 



where 



d9 r N5Sa* , N3Sa* d “a , 

dn = +^— — (e+N4)J/GW ,9) 

a a) 



G(w ,0) 

a. 



N2 (T -T ) N3 (T -T ) ,_ . 

+ _E ^£f-(— -1) (0+N4) 



Fa 



dT 



Sa 



dT w. 



+ N3 (Sa*) (— - 1) 

CO 

a 



(C) 



V* 



Equation (III. 22), V* (w 

a 



= (^u -N7(/0dn)w -u 

q ri a a 



n 



( ■ , d0 d0 
■J 9dn 'dn'dn 



d0_ 
a dn 



0 



) / (N8H* 



fg 



) 

o 



Sa*) . 



The boundary conditions for this set of equations are the 
following : 
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(a) For nondimensional temperature: 



n = 0.0 



n = 1.0 



d_0 

dn 



o 



a 




9 

c 



e = l.o 



The constant a is not known but, based on laboratory data, can 
be estimated and through an iteration process the correct value 
can be obtained. A proposed iteration scheme is presented in 
Chapter IV. 

(b) For working fluid vapor mass fraction: 

n = 0.0 w = to 

a a„ 

0 

n = l.o to =1.0 

a 




The constant 3 is not known, but can also be estimated for an 
initial value and, through an iteration scheme, a correct value 
can be obtained. 

(c) For working fluid vapor dimensionless velocity: 

n = o.o v* = o.o 
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E. AUXILIARY EQUATIONS 



A means to derive saturated working fluid heat of vapori- 
zation, density, and pressure is required before a numerical 
solution is attempted. Using methanol as the working fluid and 
data from Eubank [Ref. 8] for saturation conditions, the 
following functions are derived. 

1 . Polynomial for the Working Fluid Heat of Vaporization 
Tabulated values of the heat of vaporization for methanol 

are plotted on Figure 3.1 with a solid line. A Fortran program, 
polyfit, is used to obtain the coefficients of the following 
polynomial : 

H (T) = aT 2 + bT + c (III. 23) 

where a = -9.96 x 10 ^ Btu/Lbm ( °R) 2 , b = .730 Btu/Lbm(°R), and 
c = 397.99 Btu/Lbm. The values of H^^(T) using this polynomial 
are plotted on Figure 3.1 with a dotted line. 

2 . Polynomial for the Working Fluid Vapor Saturation 
Density 

Tabulated values for methanol saturated density are 
plotted on Figure 3.2 with a solid line. Again the Fortran 
program, polyfit, is used with the following modification: 

m T* 

T* = - where T c = 500°R and £ = e 

T _ ref 

ref 

A polynomial in the variable £ is used to describe the 
saturation density. 
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SATURATION CH OH HFG 



o 




(wai/nia) ojh 
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Figure 3 . 1 . Graph of Methanol Heat of Vaporization Versus Temperature 



SAT. Cli OH VAPOR DENSITY 



o 




(E^u/wai) vs 
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Figure 3.2. Graph of Methanol Saturation Density Versus Nondimensional Temperature 



The polynomial, Sa(£), is then: 



Sa(£) = a? 3 + b^ 2 + cE, + d (III. 24) 

where a = .1293 Lbm/ft 3 , b = -1.0475 Lbm/ft 3 , and 
c = 2.8868 Lbm/ft 3 and d = -2.7037 Lbm/ft 3 . The value of 
dSa/dT is also required for the numerical solution. This is: 



W - < 3a ? 3 + 2b « 2 + C 5 )/ T ref 



(III . 25) 



3 . Polynomial for the Working Fluid Vapor Saturation 
Pressure 

Tabulated values for methanol saturation pressure are 
plotted on Figure 3.3 with a solid line. Polyfit is used with 
the same modification. 

The polynomial Fa(£) is: 

Fa(C) = a ? 3 + b £ 2 + c£ + d 

where a = 24.5994 Lbf/in 2 , b = -191.2472 Lbf/in 2 , 
c = 498.4940 Lbf/in 2 , and d = -439.9883 Lbf/in 2 . The first 
derivative of Fa(£) with respect to temperature is equal 
to : 



dFa 

dT 



(3a C 3 + 2b£ 2 + c£)/T ref . 



(Ill . 26) 
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SAT. CTLOH VAPOR PRESSURE 
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IV. DISCUSSION AND RECOMMENDATIONS 



A. DISCUSSION 

The final version of the model for a variable conductance 
heat pipe with a vertical orientation is a system of three 
ordinary, non-linear, coupled, differential equations. These 
three equations are (III. 14), (III. 19), and (III. 22). 

The remaining challenge for the numerical model is to pro- 
pose an iteration scheme that can solve this system of three equa- 
tions. A special thank you goes to the collaboration of Dr. M.D. 
Kelleher and Dr. K.T. Yang for the following iteration technique. 
The steps of the iteration procedure are as follows: 

(a) To start the iteration process, assume a form for the 

relationship between 6 and to such that to = C-, 9 + C 9 . This 

a a -L z 

linear relationship is reasonable because: 



(i) When g = 0, both to and 9 are values close to 

a 



zero . 



(ii) When g = L both to and 9 are equal to 1.0. 

a 

(iii) This step is required for the first iteration 



only . 



The equations for this step are then: 



a 



C-^ 9 + C 2 



(IV. 1) 



dto 




a 



(IV. 2) 



dg 
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Then, using Equation (III. 19), solve for This establishes 

an 

a first degree equation, in terms of 0 only which can then 

be integrated numerically to establish a first iteration 
relationship for 0 (n ) . 

dfl 

(b) With this relationship for 0 and defined from step 
(a) above. Equation (III. 22) can be solved for V*. 

(c) Once step (b) is completed, Equation (III. 14) can be 
solved for dco /dn. The integral of dco /dg is then used to 

cL cl 

define . The process can then be repeated. 

In summary, the first iteration goes through steps a-c 
and the remaining iterations go through steps b-c until there 
are no further changes in the form of 0, w , and V*. In 
theory, the outcome of this numerical iteration is a table 
which gives dimensionless temperature, 0; velocity, V*; and 
mass fraction of the working fluid vapor, w , versus non- 

3 . 

dimensional condenser length, g. 



B. RECOMMENDATIONS 

1. This numerical analysis was written to yield results via 
a computer program. The Continuous System Modeling 
Program III [Ref. 9] is a computer package that could 

be used to solve the system of equations. The CSMP 
[Ref. 9] includes a powerful integral statement which 
carries out a Fourth Order Runge-Kutta Method. 

2. This model assumes constant values for the binary mass 
diffusion coefficient and the convection heat transfer 
coefficient. Once numerical values are obtained from 
the iteration scheme, the accuracy can be improved by 
using analytical expressions instead of constants for 
the binary mass diffusion coefficient [Ref. 10] and the 
convection heat transfer coefficient [Ref. 11]. 

3. An attempt could be made to model the variable conduc- 
tance heat pipe that is not in a vertical orientation. 
Kelleher [Ref. 3] shows that a stratification of the 
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working fluid and noncondensible gas occurs when the 
heat pipe is in the horizontal position. This position 
makes the transport processes inside the pipe three 
dimensional, which would require a finite difference 
computer model for solution. 
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